Orientational order in two dimensions from competing interactions at different scales 
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We discuss orientational order in two dimensions in the context of systems with competing isotropic interac- 
tions at different scales. We consider an extension of the Brazovskii model for stripe phases including explicitly 
quartic terms with nematic symmetry in the energy. We show that leading fluctuations of the mean field nematic 
solution drive the isotropic-nematic transition into the Kosterlitz-Thouless universality class, i.e. these systems 
have a thermodynamic phase with orientational quasi-long-range order. 
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I. INTRODUCTION 



Systems with competing interactions at different length 
scales are common in nature'. Examples go from highly 
correlated quantum systems like quantum Hall samples^ and 
high Tc superconductors^, to classical systems like ferromag- 
netic ultrathin films^, diblock copolymers^, colloidal sus- 
pensions^, ferromagnetic garnet films'' and magnetic fluids^^i. 
The essential phenomenology of these kind of systems was 
described in a classic paper by Brazovskii". Competition 
on different scales gives rise to ordered phases dominated 
by a non-zero wave vector in reciprocal space, as opposed 
to the usual k = Q long range order. The non-zero value 
ko of the dominant wave vector gives rise to spatially mod- 
ulated structures. In three dimensions Brazovskii showed that 
striped phases appear through a first order phase transition in- 
duced by fluctuations". Stripe patterns show both positional 
(anisotropic) and orientational long range order, although the 
stripe solutions in the self consistent Hartree approximation 
are marginally stable in three dimensions. Subsequent works 
extended the original results to two dimensional systems in 
spite of the fact that, strictly speaking, fluctuations prevent 
any long range order'^''-*. 

Motivated by recent experimental observations of phases 
with modulated order in two dimensional systems, we ana- 
lyzed in a recent letter'** the conditions for the existence of 
a purely orientational phase, a nematic phase, in models of 
the Brazovskii class. In the framework of the Renormaliza- 
tion Group, we showed that an isotropic-nematic phase tran- 
sition is generically present in these kind of models, pro- 
vided suitable quartic interactions between the basic degrees 
of freedom are taken into account. These interactions are 
naturally generated in the renormalization process. Further- 
more, we found that, in two dimensions, renormalization of 
the Brazovskii model gives rise to an infinite number of rel- 
evant terms which makes the model non-renormalizable. We 
have shown in Ref. 14 that all those terms possess a com- 
mon symmetry under rotations by tt, a nematic symmetry. 
Keeping only the term with the highest symmetry, correspond- 
ing to quadrupole-quadrupole interactions, we showed that a 
isotropic-nematic phase transition is present and that it is of 
second order at mean field level. Nevertheless, it was antic- 



ipated that the nature of the transition would probably be af- 
fected upon inclusion of fluctuations''*, since it is not possible 
to break a continuous symmetry in two dimensions with short 
ranged interactions'^. 

The original Brazovskii model in three dimensions is at its 
lower critical dimension, with fluctuations in the stripe solu- 
tions diverging logarithmically with the linear size of the sys- 
tem. The situation is more delicate in two dimensions where 
fluctuations are linearly divergent'^. Then, if some order of 
this kind survives in two dimensions it must be purely ori- 
entational. However, it is still necessary to check whether 
orientational order survives to fluctuations of the relevant or- 
der parameter. Of course, in real systems, the nematic or 
even the stripes phase can be stabilized by other factors, Uke 
anisotropics coming from the lattice substrate, the presence of 
impurities or some disorder that pin the stripe order. These 
effects will not be considered here, where we rest at the level 
of a completely isotropic system. 

In the present work we pursue the analysis of a general- 
ized Brazovskii model which takes into account quadrupolar 
interactions in two dimensions. We briefly review the mean 
field treatment of Ref. 14, and we evaluate thermal fluctua- 
tions. We show that the isotropic-nematic transition belongs 
to the well known Kosterlitz- Thouless universality classic, 
i.e. upon inclusion of order parameter fluctuations, the mean 
field solution with nematic long range order in fact retains 
only quasi-long-range orientational order. In turn, the stripe 
solution is unstable to fluctuations and a possible smectic-like 
phase reduces to a point at zero temperature. Similar results 
were found long time ago by Toner and Nelson in the con- 
text of defect mediated melting in two dimensions '^. In fact, 
both approaches are complementary and consistently lead to 
the same phase diagram. The present results were briefly an- 
ticipated by us in a reply '^ to a comment'^ to Ref. 14. 

In the following, we introduce in ^ a prototypical model 
of the Landau-Ginzburg type for a nematic order parame- 
ter. We show that, while the mean field treatment leads to 
a second order phase transition in two dimensions, low en- 
ergy fluctuations diverge logarithmically, as in the XY model 
of magnetism, destroying the long range order and leading to 
an algebraic decay of the correlations. In jllll we introduce 
the extended Brazovskii model considered in Ref. Q, and 
briefly discuss the mean field solution. We show that the rele- 
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vant leading order fluctuations can be mapped to an XY type 
model. We thus show that, within the Gaussian approxima- 
tion, the isotropic-nematic phase transition is in the Kosterlitz- 
Thouless universality class. The analysis also allows us to ex- 
press the effective elastic constant K{T) as a function of the 
parameters of the original model. A brief discussion of the 
results is given in ^IVI 



energy for a constant tensor near the transition, where < Q > 
is very small. The leading rotational invariant terms are: 



1 



1 



FiQ) = -a2TriQ') + -a^TriQ^) + 



(6) 



Using Eq. (|3]l it is very simple to show that the free energy 
reduces to: 



11. LANDAU-GINSBURG THEORY FOR THE NEMATIC 
TRANSITION 

The Landau-Ginzburg theory for a three dimensional ne- 
matic phase is well known^". Here, we briefly restate it for 
a two dimensional system because it presents some charac- 
teristics exclusive of the dimensionality of the problem. This 
analysis will be also helpful as a guide to the evaluation of 
fluctuations of the specific model of section HIH 



A. Order parameter 

The order parameter for a 2d nematic has two components 
to identify an orientation in the plane and the intensity, but not 
a direction. This means that it must be symmetric under the 
transformation 6 ^ 6 + tt. 

Consider a complex number written in the form: 



Q 



a e 



120 



(1) 



where 6 is an angle in 2d space. Then, if < Q 0, the 
phase is said to have orientational order in the 9 direction with 
the nematic symmetry 9 ^ 9 + tt. The nematic symmetry 
implies that the order parameter is not a vector Instead, we 
can arrange the real and imaginary parts of ([T]i in a second 
rank symmetric and traceless tensor in the the following way: 



TZe{Q) = a cos(26i) = Q^x = -Qyy 
Im{Q) = a sin(26') = Qxy = Qyx 



(2) 



and 



^ j Q XX Qxy 

Qxy Q XX 



cos(26') sin(20) . 
" ' sin(26i) - cos(26l) ' 



Defining now a unit vector (the director) fi with components 
n^, — cos 9 and Uy = sin 9, the nematic order parameter 
reads: 



F(a) 



1 



1 



-a2a 



(7) 



Note that at this level the free energy is independent of 9, that 
means that it is invariant under arbitrary global rotations. In 
particular, in two dimensions, Tr{Q^) — and therefore there 
are no terms with a^, at variance with the id case. This im- 
plies that the mean field isotropic-nematic transition is of sec- 
ond order in 2d. 



B. Mean field phase transition 

Consider the free energy of Eq. (|7]i and suppose that 04 > 
0. Therefore, if 02 > the only minimum of this energy is 
a ~ Q and then < Q >= 0. Conversely, if 02 < the 
minimum is at a = (—02/04)^/^ and 



< Q >= 



""2 A2S 
04 



or in terms of the director components: 



< Q., 2 



-Q2 
a4 



2 ^ 



(8) 



(9) 



To leading order (near the transition) 02 (T) — a {T — T*), 
where a > is a constant and T* is the critical temperature. 
At the critical point the rotational symmetry in the plane is 
spontaneously broken. Choosing the director direction to cor- 
respond to 6* = then: 







if T > T* 



(Q) 



I a ^J^* 2^^1/2 jf y <^ J^* 



(10) 



This is the classic Landau-Ginzburg scenario for a second or- 
der phase transition with Tc — T*. We will see that fluctua- 
tions in the director orientation change this picture. 



a 



nl - riy 2nxny 



2nxny Uy - 



or in component notation: 



Qij = 2a ( rii - -n'^dij 



(4) 



(5) 



Equations ([T]i and Q are two different ways of writing the 
same thing. We can now developed a Landau-Ginzburg free 



C. Fluctuations 

We have developed the free energy of Eq.© considering 
that the order parameter Q is constant. However, if we want 
to study local fluctuations we can consider a local order pa- 
rameter of the form Q = Q{x) and study the free energy for 
small variations of Q{x) around the mean field value Q. In 
order to do this we need to introduce terms proportional to 
derivatives of the order parameter in the expansion of the free 
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energy . To leading order, we consider just first derivatives of 
Q and write a rotational invariant free energy of the form: 

F{Q) ^7/^'^ {^TriQDQ) 

+ Ja2Tr(Q2) + ia4rr(Q4) + ...| (H) 

where p is a stiffness constant and the symmetric derivative 
tensor Z)y = ViVj. 

Because the free energy is symmetric under global rota- 
tions, low energy angle fluctuations are the most relevant 
modes that rule the behaviour of the system. We will see that 
in 2d the angular correlations are logarithmically divergent, 
ruling out true long range order but showing instead quasi- 
long-range order or power-law decay of spatial correlations. 
Consider a local order parameter of the form 

V 04 

or in tensor form, as a function of the director components: 



Qij{x) = 2^—^ yni{x)nj{x) ~ -rfiSijj (13) 

Thus, fixing the modulus to its mean field value, we proceed 
to study small local fluctuations in the direction of the director 
in the nematic phase. Replacing Eq.(fT2b or ( fT3] l into ( fTTT i we 
find that: 

5F = F{Q{x)) - F{< Q>)= K{T) j d?x \Ve{x)\'', 

(14) 

where 

K[T) (15) 

Therefore, the free energy for the small angle fluctuations of 
the director can be mapped into the free energy of the XY 
model Angle correlations in the XY model decay alge- 
braically as: 

(cos(6l(x)-e'(0))) cxa;^" (16) 

with 77 = T/2TTp. Then, the isotropic-nematic transition in 
2d belongs to the Kosterlitz-Thouless universality class with 
a disordering mechanism mediated by the unbinding of topo- 
logical defects The only difference is that the role of vor- 
tices in the XY model is played here by disclinations'^. 

These results are independent of any microscopic mecha- 
nism. In the present case, if we begin with a Brazovskii type 
Hamiltonian, one should be able to reach Eq. (fTTl i where the 
parameters 02, 04, p and T* should be written in terms of the 
more "microscopic" onesii. This is the subject of the next 
sections. 



III. MODEL WITH COMPETING ISOTROPIC 
INTERACTIONS AT DIFFERENT SCALES 

Long time ago Brazovskii" introduced a rather general 
model with the aim of capturing the physics of systems 
with isotropic competing interactions at different scales. The 
model should be relevant for a wide class of systems as dis- 
cussed in the Introduction. Specializing to two spatial dimen- 
sions and considering a scalar order parameter (Ising symme- 
try ), the Brazovskii model is defined (in reciprocal space) by 
a coarse-grained Hamiltonian of the type: 

(17) 

where ro(r) - a{T - T^), k = |fc| and /cq ^ \ko\ is a 
constant given by the nature of the competing interactions. 
J^d'^k = J^^ dO S^'^^^ dk k and A ~ ^/mro is a cut-off 
where the expansion of the free energy up to quadratic order in 
the momentum makes sense. The "mass" m measures the cur- 
vature of the dispersion relation around the minimum fco and 
the ellipsis in eq.(fT7]i indicates higher order terms in (fc — fep)- 
The correlator has a maximum at fc = fco with a correlation 
length ^ ^ 1/ ^/mfo- Therefore, near criticality (rg — > 0), 
the physics is dominated by an annulus in momentum space 
with momenta fc ~ fco and width 2A. This implies that at 
high temperatures the model possess a continuous symmetry 
in momentum space, or in other words, a large phase space 
for fluctuations. The original model proposed by Brazovskii 
contains also an interaction term proportional to 0"*. In the 
mean field approximation, this model leads to a second order 
phase transition from an isotropic phase at high temperatures 
to an anisotropic stripe phase with modulation of the order 
parameter in the form: 

{(l){x)) = A cos{kox). (18) 

Working in three dimensions, Brazovskii showed that includ- 
ing fluctuations of the order parameter self-consistently leads 
to a "fluctuation-induced first order transition". Subsequently 
this transition was observed and studied in diblock copoly- 
mers^' . In two dimensions, stripe phases arising from compet- 
ing interactions are also observed in many systems, a notable 
example being that of 2d ultrathin ferromagnetic films with 
perpendicular anisotropy, in which the short range exchange 
interaction between spins is frustrated by the long range char- 
acter of the dipolar interaction, giving rise to the well known 
magnetic domains^^. In recent years there have been indica- 
tions that a mechanism similar to that proposed by Brazovskii 
can be at work in these systems^^. However, stripe solutions 
are not stable with rigorously isotropic interactions, and fluc- 
tuations in the stripe direction diverge logarithmically in 3d. 
Then, for two dimensional systems the situation should be 
worse unless some isotropy-breaking effect be at work, like, 
e.g. lattice effects'^. Nevertheless, even if positional long 
range order is forbidden for such models in 2d, one can ask if 
some kind of orientational order, reminiscent of stripe order, 
may survive in an isotropic model of the kind considered. 
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Recently'"*, we analyzed which kind of interaction terms 
could give rise to purely orientational order in two dimensions, 
besides the already known Brazovskii stripe solutions, which 



posses orientational as well as translational long range order. 
Considering a generic interaction term of the form: 



^^int — 



(Pki (fk2 (Ph d^h 
A (2^(2^(27r)2 (27r)^ 



u{ki,k2, fcs, (t>{ki)<j){k2)<j){k3)(j}{k4) 5{ki + k2 + ky, + k^). 



(19) 



and performing an analysis in the context of the Renormal- 
ization Group, we found^i that as a consequence of the 
isotropy of interactions in 2d the space dependent function 
u{ki, k2,k3,k4) should depend solely on an angle 6 and have 
the form: 



u{6) = uo + U2 cos(26') + U4 cos(46') 



(20) 



Here, the first term uq leads to the usual 0"' theory considered 
by Brazovskii in his model for the isotropic-stripe transition. 
The other terms are all relevant in the RG sense and then the 
theory is not renormalizable. Nevertheless, it is evident that 
all the terms share the nematic symmetry 6 ^ 9 + tt. We then 
proceeded to analyze the effect of the first of those terms, pro- 
portional to cos (29). From the definitions in section UTaI we 
realize that the cos (26*) factor can be conveniently expressed 
in terms of the tensor order parameter and the interaction en- 
ergy can be written in the form: 



Hint — J <fx ^uq (f)"^ {x) + U2 tr + tr 
with 7 > and 



Q^} (21) 



1 



^) V,V, - (/.(x). (22) 



The gradients are related to the director n,; — Vi/|Vi|. From 
Eq.(l22]i it is clear that the nematic order parameter is essen- 
tially a quadrupolar moment. 

Next, we proceed to analyze this extended model in the self 
consistent Hartree approximation. 



A. Hartree approximation 

Replacing in eq. (1211 1 cf)'^ cjp' < (f)^ > and tr 
tr {0 (VjVj — ^V^f^y) (/)} {Qij{x)), where the mean val- 
ues have to be determined self consistently, we obtain a 
quadratic Hamiltonian in the Hartree approximation, which 
in reciprocal space reads: 



Haitic 



(Pk 

(2^ 



m (p-^C-\k)) ^{-k), (23) 



with the two-point correlator C{k) given byii: 

T 

C{k) 



r + i(fc - koY - a^k^ cos(26')(u2 + ja^) 



(24) 



Here 



and 



r = To + Uq 



_fk_ 

(2^ 



C(fc) 



_fk_ 

(2^)2 



k^cos{2d) C{k) 



(25) 



(26) 



' is the 



where we have chosen (Qij) — 2a {rii rij — ^nl^Si 

angle subtended by k with the director h. 

Equations ( l24l i. dZST l and ( l26b must be solved self- 
consistently. Its solution has been discussed in Ref. [3 The 
main result which comes out is that in the case of attrac- 
tive quadrupole interactions, U2 < 0, equation ( |26l ) has non 
trivial solutions for the nematic order parameter. Writing 
the equations in terms of adimensional parameters r rT, 
ko fco \/ mT and tq t = a{l — Tc/T), one finds out that 
for high temperatures, T > Tc, the only possible solution is 
a = 0. Nevertheless, atT — Tea nematic phase emerges con- 
tinuously with a ~ c(Tc — TY^^ and the critical temperature 

A spontaneously broken continuous symmetry in a two di- 
mensional system with short range isotropic interactions is 
forbidden by the Mermin- Wagner theorem'^. In these sys- 
tems, fluctuations of the order parameter typically diverge and 
the precise nature of the divergence can say if order is lost 
exponentially fast or if it decays more slowly giving rise to 
what is called "quasi-long-range-order". In order to analyze 
the effect of fluctuations of the nematic order parameter near 
the transition, we write the free energy of the model in the 
Hartree approximation. The partition function is 



V(t> e 



-I3Fh 



Integrating over one arrives at 



1 



(27) 



(28) 



In the hmit [k — k^Y /m << Vc (that is very near the Hartree 
critical temperature), the free energy reads 



1 T 1 f^f jk-ko? 



a^A:^U2 cos(26l) 



(29) 
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where Tc and a are the solutions of the self-consistent Hartree 
equations ( |25] ) and ( |26] l described in Ref . [l^. 

Since near the transition a^kQU2/rc << 1, we can expand 
the free energy in the form: 



H 



'''"'-Tr{fc2cos(20)} - ^Tr{fc4cos2(20)} , 



(30) 



where a fc-independent term was absorbed in Fh- The first 
term of the expansion is zero by symmetry (upon integration 
over 9), therefore at leading order the free energy reduces to 



Fh 



4/3r 



-TT{k 



4cos2(2(0-(p))} 



(31) 



This expression represents the contribution of the anisotropic 
(nematic) part in the Hartree approximation, very near the 
transition into the nematic phase (a ^ 0). In the last expres- 
sion we have introduced the angle that is the reference from 
which we measure the angle 6 (which is the integration vari- 
able). At this level is an arbitrary constant, as it should be 
in a spontaneous symmetry breaking scenario. Next we will 
study smooth fluctuation of this field. 



B. Fluctuations 

In the same way we have done in the Landau-Ginzburg the- 
ory (and for the same reasons), we consider angle fluctuations 
of the order parameter, Q{x) ^ a exp i2iy9(a;). Therefore, the 
free energy now takes the form: 



Fh = - 



4:f3r; 



:TT{k'^ cos^ 2{0 ~ (p{x))} 



(32) 



The difficulty with this expression is the evaluation of the 
trace, since its argument is not diagonal neither in k nor in 
X space. To evaluate it, we make a coarse graining of config- 
uration space, in such a way that in a small region around a 
point xq we consider ip essentially constant. Then, we can av- 
erage over all points xa covering all configuration space. This 
coarse grained free energy can be diagonalized in k space and 
the trace can be easily evaluated. Consider the following ex- 
pansion for ip{x), for a fixed point xq: 

ip{x) = (p{xq) + V ip{xo) ■ {x - Xq) + . . . 

« ip'{xo) + \'ip{xo) ■ X (33) 

where the constant (p'{xo) = fixo) — W(p{xo) ■ xq. With this 
expansion we rewrite the cosine in the expression for the free 
energy: 

cos 2(0 — cos 2(0 — </3q — V(^(a;o) • x) 

« cos2e' + 2{Vtp{xo) ■ x) sin 261' (34) 



where 6' = 6 — (pg and we have considered smooth fluctua- 
tions |V(^(a;o)| << 1. Therefore, 

cos^ 2(6* - (/?) w cos^26l' + 4cos26l'sin26'' Vv3(a;o)-x 



+ 4(V(y9(xo) • xf sin^ 29' 



(35) 



The first term contributes with an additive constant to the free 
energy, and then we will not consider it anymore. The sec- 
ond term is identically zero by symmetry considerations, as 
shown in the appendix. The relevant leading contribution to 
the fluctuations is the last one. Thus, let us consider the coarse 
grained free energy for smooth fluctuations: 



F: 



fl 



-4-9 f i9 
a U2 / a Xq 



V 



Tr 



|fc4(V(^(xo) • xf sin2 26''| , 



(36) 

where V is the volume of the system. Using the representation 
X — iV k, we write the trace in k space in the form: 

Q,4y2 r r dk ^ ^ 

Pfi = J d^^o J -^kd9 sm^29iVip{xoyVk fk\ 

(37) 

The k derivatives can be evaluated as 

(yp>{xa) ■ Vkfk'^ = 4\Vipo\^k^ {l + 2sh? 9] , (38) 

where V(^o • k = |V(^o||fc| sin0 because Vpo is in the direc- 
tion of the fluctuations (perpendicular to the director) and then 
the angle between V(^o and k isTr/2 — 9. 
We finally obtain 



(39) 



where 



r = 4 / fc3 / gij^2(20) {l + 2sin2 0} 

Ao-A (27r)^ Jo 



TT 



(40) 



to leading order in the cut-off A ~ y/mro = ^Jma{Tc — T). 

Remembering that the Hartree solution of the order param- 
eter is a = c{Tc — T)^/^, we can write the free energy for the 
fluctuations as: 



Fn = K{T) J d^x \Vip{x)\^ 
where the elastic constant is given by: 



5/2 



(41) 



(42) 



The constant k = {8/TT){a^/^c%^\l^'^)/{my^ul^^^). Ex- 
pression ( l4n i is equal to equation (fT4l i. Then, Gaussian fluc- 
tuations of the nematic order parameter around the mean field 
solution diverge logarithmically, and the nematic phase does 
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not have true long range order, but instead retains quasi-long- 
range order with the well known Kosterlitz-Thouless phe- 
nomenology"'. The conclusion is that fluctuations change 
the nature of the phase transition and, in particular, the crit- 
ical temperature departs from its mean field value Tc- The 
isotropic-nematic phase transition in the present model takes 
place at a temperature Tkt- At this temperature, a continuous 
phase transition mediated by unbinding of disclinations hap- 
pens with Tkt = iT^/S>)K{TKT)^- This relation, together 
with (l42l i. allows to obtain the transition temperature as: 



Tkt = 



1 



1 + ^ 

7TK , 



(43) 



IV. CONCLUSIONS 

Systems with competing interactions can develop complex 
ordered phases , with characteristics different from the usual 
ferromagnetic long range order For many systems, competi- 
tion may lead to ground states with modulations in the order 
parameter These broken symmetry states naturally show ori- 
entational order, and sometimes also positional one. While 
positional long range order is strongly suppressed in two di- 
mensions for systems with isotropic interactions and contin- 
uous symmetry, orientational order is more robust. We have 
studied a rather general model for competing interactions at 
different scales, looking for conditions for the existence of a 
purely orientational phase at low temperatures. 

We have shown that, in the two dimensional Brazovskii 
model, the quartic interactions with higher derivatives of the 
order parameter are all relevant terms in the Renormaliza- 
tion Group sense. All these terms can be arranged and inter- 
preted as representing multipole interactions. Among them, 
the quadrupole-quadrupole interaction is the first non trivial 
contribution. A mean field solution of the model gives rise an 
isotropic-nematic phase transition. The analysis of Gaussian 
fluctuations around the mean field solution leads to a phase 
diagram similar to the one found by Toner and Nelson in the 
context of defect mediated melting in two dimensions—. Po- 
sitional order of Brazovskii stripe solutions is destroyed by 
thermal fluctuations, which are known to diverge linearly in 
two dimensions. However, orientational quasi-long ranged or- 
der is preserved in the nematic phase of the extended model. 
We have shown that there is a critical temperature Tkt at 
which orientational quasi-long-range order is destroyed. By 
analogy with the XY model, one can conclude that the disor- 
dering of the nematic phase takes place by means of a discli- 
nation unbinding mechanism, and the isotropic-nematic phase 
transition is in the Kosterlitz-Thouless universahty class. 



The main difference of our approach with that of the Toner- 
Nelson-Kosterlitz-Thouless is that our model allows for an 
analysis of both sides of the phase transition. This fact makes 
possible to characterize the transition in terms of "micro- 
scopic" parameters, which describe the underlying competing 
interactions. Also, within the present formahsm, it is possible 
to alternatively interpret the nematic phase as a quadrupole 
condensation rather than a melting of topological defects. 

Finally, the presence of a nematic phase from competing 
interactions in two dimensions can be present in a variety 
of systems like ultrathin ferromagnetic films with perpendic- 
ular anisotropy^"*, block copolymers'^, microemulsions and 
colloids*''^'', between others. The detection and quantitative 
characterization of such phases in those systems rely on novel 
imaging techniques which are at present rapidly evolving. 



V. APPENDIX 

We show here that the second term in eq. ( |35] ) gives zero 
contribution to the free energy. The contribution of this term 
to the free energy is: 



F2 ~ 4i / d^fccos26''sin26'' \I^}{xq) ■ Vfefc^ 



(44) 



with 



V(p(xo) • Vfcfc'' = 4fc2v^(a;o) • fc = fc^|V(/3(xo)| sin6'. (45) 

The last expression is due to the fact that if we measure Q 
from the director and consider that V(y3(xo) is orthogonal to 
it, then the angle between Vi^(a;o) ™d fc is 7r/2 — Q. Then 
the scalar product is written in terms of cos(7r/2 — Q)= sin Q. 
Introducing in Eq. (|44] | and remembering that 9' = 9 — (p^, 

F2 - 16i|V^(a;o)| x 

X J dkk^ d9 cos2{e — ip'q) sin 2{9 — ip'q) sin 9. 

(46) 

The angular integral is identically zero whatever the value of 
ipQ, thus F2 = . 
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